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We have measured individual tunneling events and Coulomb step shapes in single-electron boxes 
with opaque superconductor-normal metal tunnel junctions. We observe anomalous broadening of 
the Coulomb step with decreasing temperature in a manner that is consistent with activation of first- 
order tunneling events by an external dissipative electromagnetic environment. We demonstrate that 
the rates for energetically unfavourable tunneling events saturate to finite values at low temperatures, 
and that the saturation level can be suppressed by more than an order of magnitude by a capacitive 
shunt near the device. The findings are important in assessing the performance limits of any single- 
electronic device. In particular, master equation based simulations show that the electromagnetic 
environment realized in the capacitively shunted devices allows for a metrologically accurate charge 
pump based on hybrid tunnnel junctions. 



I. INTRODUCTION 



Various kinds of electron pumps and turnstiles based 
on single electronics are under active investigation at 
the moment in order to redefine the unit ampere and to 
provide the current source for the quantum metrological 
triangle experimenti. The performance requirements in 
metrology are demanding: To improve the existing un- 
certainties of fundamental constants, elementary charges 
have to be transferred at a frequency of the order of 
1 GHz with an error rate less than 10"^. With present- 
day thin-film technology, it is straightforward to realize 
electron pumps where the error rates due to thermally ac- 
tivated one-electron processes are below the metrological 
bound at sub-Kelvin temperatures. Residual pumping 
errors in normal metal pumps are generally attributed 
to co-tunneling^ and environmental activation (EA)^, al- 
though a quantitative theoretical analysis of the exper- 
imentally observed error rates is still lacking. A more 
recently developed candidate is the hybrid turnstile, con- 
sisting of a single-electron transistor (SET) with normal 
metal-insulator-superconductor (NIS) junctions^. Theo- 
retically, the SINIS turnstile has the potential to reach 
metrological accuracy^, but the experimental results have 
been limited by the leakage current present when the 
turnstile is biased at its operating point^. Leakage in NIS 
structures is often attributed to a finite density of states 
(DOS) within the Bardeen-Cooper-Schrieffer (BCS) en- 
ergy gap, which is modeled with the phenomenological 
Dynes DOS>^. In Ref. d, it was shown algebraically that 
a high temperature environment (fceTcnv ^ A) can man- 
ifest itself as smearing of the DOS in the superconducting 
electrode of a bare NIS junction, which was experimen- 
tally verified. On the contrary, in a Coulomb blockaded 
system embedded in an environment, the device response 
depends on one-electron tunneling rates in a way that 
cannot be expressed in terms of an effective DOS. A tun- 
neling electron can overcome the Coulomb blockade by 
absorbing energy from the environment, which degrades 
the performance of the device. 



In the present work, we study experimentally hybrid 
tunnel junctions in the EA regime. We demonstrate that 
EA causes the observed low-temperature saturation of 
one-electron tunneling rates near the degeneracy point 
between two charge states, in the energy range that deter- 
mines the accuracy of the SINIS turnstile when operated 
with a sinusoidal drive signal. Hybrid tunnel junctions 
are well-suited for this investigation, as the EA effects 
are dramatic owing to sharp features in the supercon- 
ductor DOS, which has been previously exploited in the 
study of photon-assisted tunneling by microwave irradi- 
ation in SETs^. The simplest system where a quantized 
island charge can be observed is the single-electron box, 
first studied in Ref. 10- In the present work, we probed 
the charge state of a galvanically isolated hybrid single- 
electron box by a capacitively coupled SET electrometer. 
Owing to the superconducting energy gap, the tunneling 
rates near degeneracy are sufficiently low so that individ- 
ual tunneling events can be observed without additional 
trapping nodes that are required for normal metal de- 
vices. We emphasize that this is the same qualitative 
difference that allows quantized current pumping using 
the SINIS turnstile with only one Coulomb-blockaded is- 
land, whereas the normal metal pumps always have at 
least two islands. For demonstrating the EA effects, we 
measured samples with and without a ground plane, i. e., 
a capacitive shunt in the electrical leads acting as a filter 
for high-frequency environmental noise. We present data 
from four samples in total: Samples Gl and G2 have a 
ground plane as shown in Fig. [ija), whereas reference 
samples Rl and R2 do not have one. Samples Gl, G2, 
and Rl were fabricated simultaneously on the same chip 
and have identical design for the metal layers. Sample 
R2 is otherwise similar, but has a SISIS type detector as 
opposed to the SINIS type found in all the other samples. 
Figurelllb) illustrates the layout of the active region with 
the single-electron box and the electrometer. 

Two gate electrodes allow for independent tuning of 
the detector and the box. Whenever the box gate posi- 
tion was changed, a compensation signal was subtracted 
from the detector gate voltage canceling the direct capac- 
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FIG. 1. (a) Large scale electron micrograph of a sample show- 
ing the bonding pads and the electrical leads to the active re- 
gion in the middle. A copper ground plane isolated from the 
thin-film metal strucutres by 120 nm aluminium oxide layer 
is indicated in the bottom half of the picture, (b) Active re- 
gion of a single-electron box similar to the measured samples. 
The hybrid single-electron box is shown in the middle, and 
the electrometer SET below it. The thin- film copper layer 
appears lighter compared to the aluminium layer. Tips of 
the gate electrodes can be seen at the bottom and top-right, 
(c) Four-second real-time electrometer trace (black) recorded 
from sample G2 demonstrating the ability to accurately con- 
trol the charge state of the box. The control sequence (red) 
applied to the box gate is as follows: 0-1 s Degeneracy point, 
also detailed in panel (d), showing stochastic switching be- 
tween the two charge states with nearly equal occupation. 
1-3 s Sinusoidal drive between the charge states at 20 Hz, 
also detailed in panel (e). 3-3.5 s and 3.5-4 * "Hold" mode 
in the first and second charge state, respectively. The sample 
stage temperature is 60 mK. 



itive coupling to the detector island. For the following 
discussion, one can consider the compensated box gate 
voltage to be the only external control signal. Due to 
the compensation, charge states of the box correspond 
to fixed detector current levels. The response of the 
box state to a gate control sequence is demonstrated in 
Fig.[ljc)-(e), displaying a real-time measurement of sam- 
ple G2. 

All the tunnel junctions were oxidized strongly, result- 
ing in specific resistances of 3 kft /im^ , or about 300 kft 
for the detector junctions and up to 3.5 Mfl for the 
smaller box junction. High resistance is required for the 
box junction to bring the tunneling rates down into the 
measurement bandwidth. The detector current was read 
out in DC-SET configuration using a room-temperature 



current amplifier, limiting the usable bandwidth to about 
5 kHz. 



II. THEORY OF ENVIRONMENTALLY 
ACTIVATED TUNNELING 

Due to small tunnel junctions and the lack of galvanic 
connection, the charging energy of the box for a single- 
electron tunneling was at least 4 K in all of the samples. 
Hence, in all of the experiments performed in dilution 
refrigerator temperatures, at most two charge states of 
the box have a non-negligible occupancy. We can safely 
assume that an equilibrium quasiparticle energy distribu- 
tion is reached before each tunneling event as the highest 
tunneling rates considered in the experiment were below 
10 kHz. Thus, a complete model for the internal dy- 
namics of the box is a two-state fluctuator, the states of 
which we will denote by 1 and 2. Neglecting back-action 
from the detector, the transition rates are symmetric and 
given by r='= = T{±E), where E = U2 — Ui is the differ- 
ence in the electrostatic energy of the two charge states, 
and r(£') is the Golden rule tunneling rateii 

T{E) = ^— / dE, / dE2ns{E,)fs{Ei) 

C JlT J —00 J —oc 

[I- h{E2)]P{Ei-E2+E). (1) 

Above, ns is the DOS for the superconducting elec- 
trode, /n(s) is the occupation factor for quasiparticle en- 
ergy levels in the N(S) electrode, and P{E') gives the 
probability for emitting energy E' to the environment. 
In our model, we use the pure BCS density of states 
ns{E) = |Re ^^f_^^ with A = 200 ^eV for the su- 
perconducting gap parameter of aluminium. As we will 
show in Appendix |B1 the quasiparticles in the N elec- 
trode are thermalized by electron-phonon coupling, and 
in the S electrode through the large overlap with the 
copper shadow. Hence, the occupation factors are taken 
to be Fermi functions at sample stage temperature Tq, 
i.e., fN,s{E) ^ [H-exp(£://cBTo)]" The probability 
for a two-electron process to change the box charge state 
is low, as for all those processes the energy cost due to 
charging energy and Cooper pair breaking is the same or 
greater, but the rate prefactor is much smaller than for 
one-electron tunneling. 



III. MODEL OF THE ELECTRICAL 
ENVIRONMENT 

Electromagnetic environment appears in the P{E) the- 
ory through the expression Re [Zt(w)], where Zt{uj) is the 
total impedance of the environment including the paral- 
lel junction capacitance, as seen from the two electrodes 
of the junction. We model the external environment as 
a parallel i?C-element with the resistor at liquid Helium 
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temperature in the same manner as in Refs. HI and[l. As- 
sume that the box electrodes have a mutual capacitance 
of Cj , the main contribution of which is from the junction 
area, and that the two electrodes are coupled to the en- 
vironment through capacitances C'l and C2 , respectively. 
Then, the impedance Zt{uj) can be evaluated as 



iujCn 



iojC 



(2) 



where Co = (1/Ci + I/C2) is the series capacitance 
of the coupling capacitances Ci and C2, and A\\B 
{1/A + 1/B)~^ is used to denote parallel impedance. Un- 
der the assumption Cj ^ C, we have the algebraic result 



Re[Ztiuj)] = Re 



R/k 



l + iujiR/K){KC) 



(3) 



where k = 1 + Cj /Cq . We refer to k as the decoupling 
factor, as K = 1 corresponds to the case where a bare 
junction is directly connected to the environment. In the 
hybrid single-electron box samples, a typical box junc- 
tion has Cj — 0.2 fF based on the observed charging 
energies of 4 — 5 K, whereas the couplings are of the or- 
der Ci_2 = 0.02 fF based on the gate modulation periods 
and geometric modeling. Hence, one can expect n of the 
order of 20. It is hard to give an accurate estimate of k 
based on the conductor geometry as the environmental 
noise can couple to the box from all of the gate and bias 
lines, but on the other hand the common-mode compo- 
nent of the noise does not affect tunneling through the 
junction. We do not try to determine an explicit value for 
K, but instead use R^s ■— R/k and C^ff ■— kC directly 
as fitting parameters. When computing the theoretical 
predictions, the P{E) function for given environmental 
parameters is evaluated numerically. 



IV. MEASUREMENT RESULTS 

A basic measurable quantity of the box is the charge 
state occupation probability that depends on the ratio of 
forward and backward tunneling rates, and assumes the 
form p = r-/(r+ + r-) for state 1. For an NIN junc- 
tion under negligible environmental influence, one has the 
analytic resuh r± = giving a Fermi 

form for the occupation as a function of the energy differ- 
ence E of the charge states, pnin = [1 + exp(ii^/fcBT'o)]~^. 
Occupation factors can be determined from the Coulomb 
staircase, i.e., the detector DC current as a function of 
the control gate, by normalizing the current levels of con- 
secutive plateaus to and 1, respectively. Because at 
most two states are active, detector non-linearity with 
respect to box charge does not affect the results. The 
occupation can be determined whenever the detector has 
sufhcient sensitivity to the box charge, even when the 
tunneling rates are beyond the measurement bandwidth. 
If the tunneling rates are low enough so that the detector 



can separate the charge states, we calculate the occupa- 
tion from the weight of histogram peaks which gives im- 
proved immunity to low-frequency drifts in the detector 
signal. 
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FIG. 2. (a), (b) Grayscale lines: Coulomb staircase in the NIS 
state, i. e., charge state occupation probability near the degen- 
eracy point as a function of the energy difference for samples 
Gl and Rl at temperature points in the range 60-200 mK cor- 
responding to the step width data presented in Fig.[lja), (b). 
The curves shown here are obtained by analyzing real-time 
electrometer traces where the two states are separable. For 
sample Gl with the ground plane, the step sharpens with re- 
duced temperature, whereas for the reference sample Rl the 
opposite is observed. Colored dashed line: Staircase in the 
NIN state at the base temperature, (c)-(f) Step widths as a 
function of temperature in NIS state (black symbols) and NIN 
state (red symbols). Black triangles correspond to step widths 
determined using real-time traces. In NIN state (red squares), 
and in NIS state at elevated temperatures (black squares), the 
step widths have been determined from electrometer DC cur- 
rent averaged over the charge state fluctuations. Dashed lines 
show the predictions of P{E) theory for NIS state (black) and 
NIN state (red). 

We determine the width of a step shape obtained from 
experiment or numerical calculation by fitting a Fermi 
function to it. Numerical results indicate that for the en- 
vironmental parameters considered in the present work, 
the Coulomb steps at high temperatures in both NIN 
and NIS states have the Fermi form given above, with 
the width corresponding to the common temperature of 
the bath and the quasiparticles in both electrodes. In 
NIN state, the linear relation between step width and 
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temperature is valid down to our base temperature of 
60 mK. However, in NIS state the step width behaves 
anomalously at temperatures below 200 mK with strong 
dependence on the environmental parameters. Experi- 
mentally, the hybrid box can be brought into NIN state 
by an external magnetic field. In Fig. [2] we show the 
data from all four samples measured in the temperature 
range 60-300 mK (60-400 mK for sample R2) in both 
NIS and NIN states. The NIN data in panels (c)-(f) 
for all samples displays linear behavior down to the base 
temperature, indicating that the box quasiparticle tem- 
peratures were not elevated due to external noise heating. 
The energy difference E depends linearly on the control 
gate voltage, with the conversion factor determined by 
the coupling capacitances and the box charging energy. 
We determined the conversion factor from gate voltage 
offset into energy by requiring that the step widths in 
the linear regime correspond to energy fee To . For sam- 
ples Gl and Rl, we display theoretical predictions using 
environmental parameters that yield best fit to the data 
presented in both Figs. [2] and [H For sample R2, provi- 
sional parameters reproducing the observed step widths 
are used. The data for sample G2 can be reproduced by 
parameters that are similar to those of sample Gl. We 
note that if at least part of the environment noise is ca- 
pacitively coupled, adding shunting capacitance over the 
sample will also decrease the effective source resistance. 

In the NIS data for the unshielded samples Rl and 
R2, EA effects can be observed at temperatures below 
200 mK as broadening of the step with decreasing tem- 
perature. At low temperatures, the tunneling rates for 
small energy differences are set by EA, and consequently 
the step width is determined by the environmental pa- 
rameters. When the temperature is sufficiently high, 
thermal activation (TA) due to broadening of the quasi- 
particle energy distribution begins to govern the tunnel- 
ing rates, and a linear relation between step width and 
temperature is restored. The EA-limited step width can 
correspond to a higher temperature than where the tran- 
sition to TA regime occurs, which results in a width mini- 
mum at finite temperature. The measured step width for 
ground plane samples Gl and G2 is nearly linear for the 
whole experimental temperature range which is reflected 
in smaller RcS and larger Coff for the environment. 

We determined the tunneling rates by stepping the 
control gate voltage across the degeneracy point, iden- 
tified as the midpoint of a Coulomb staircase step. At 
each gate position, a three-second trace of the output of 
the detector current amplifier was recorded at 100 kHz 
sampling rate followed by digital low-pass filtering with a 
sharp cut-off at 5 kHz. The whole sweep was discarded if 
an abrupt change in the detector or box operating point 
due to shifts in background charges was observed. The 
tunneling events were identified using threshold detec- 
tion, with the transition thresholds set based on the his- 
togram computed for a trace. The tunneling rates are 
estimated as F^ = T^2•j: where Ti(2) is the mean ob- 
served lifetime of state 1(2). Due to the possibility of 
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FIG. 3. (a), (b) Tunneling rates as a function of energy differ- 
ence for samples Gl and Rl, respectively, at temperatures in 
the range 60-200 mK corresponding to the NIS state data pre- 
sented in Fig. [SJa), (b). Black curves are experimental data 
whereas red curves show the theoretical predictions using the 
same parameter values as in Fig. [21 For reference, we plot or- 
thodox theory rates without environment as dashed blue lines. 
The experimental rates are most reliable near the degeneracy 
point. Gounting errors due to missed events increase as one of 
the tunneling rates grows beyond the detector bandwidth, ex- 
plaining the discrepancy between experiment and theory for 
sample Rl at high temperatures, (c), (d) Individual traces 
at the base temperature near the degeneracy point for sam- 
ples Gl and Rl, respectively, demonstrating the reduction of 
tunneling rates by capacitive shunting. 



missed events, the rate estimates obtained from a trace 
are reliable only if both rates satisfy 1/to ^ F^ ^ Sdot, 
where to is the length of the trace and i?dct is the detector 
bandwidth. 

In Fig. [31[a)-(b), we present the measured tunneling 
rates at different bath temperatures in the range 60- 
200 mK for the shielded sample Gl and the reference 
sample Rl. We also show the P{E) theory predictions 
using the same environment parameters as in Fig. [21 
For comparison, the tunneling rates without EA are also 
shown. The tunneling resistance Rt of the box junc- 
tion, an otherwise unaccessible parameter, was deter- 
mined from the data at 180 mK where tunneling rates 
depend weakly on the environmental parameters, but are 
still within the detector bandwidth. The values obtained 
this way are 3.5 Mfi for sample Gl and 2.0 Mfl for Rl. In 
the tunneling rate data, the agreement is good for the two 
samples in the E range where both rates F^ = F(±i?) 
are below 1 kHz. Above that, the experimentally de- 
termined rates are smaller than the real values due to 
missed events. For both samples, the floor where low- 
temperature rates saturate to is set by the environment. 
In sample Rl, the tunneling rates recorded at 125 and 
140 mK are lower than the base temperature rates, which 
we interpret as an artifact due to reduced sensitivity of 
the detector SET at higher temperatures. As illustrated 
in Fig.[3fc)-(d), in the ground plane sample Gl the base 
temperature tunneling rate at the degeneracy point is an 
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order of magnitude smaller than in the unshielded sample 
Rl, even after accounting for the difference in tunneling 
resistances. 



V. DISCUSSION 

To relate the observed tunneling rate data to the 
pumping performance of the SINIS turnstile, we per- 
formed master equation based simulations of electron 
transfer during the pumping cycle. For the simulations, 
we assumed a high-quality turnstile with realistic param- 
eters Rt — 0.5 Mfl per junction, Ec/k-Q = 5 K, and a 
quasiparticle temperature of 50 mK for all electrodes. 
Gate-dependent first-order tunneling rates were calcu- 
lated according to Eq. (jlj using the environment param- 
eters that we extracted for the shielded sample Gl. At 
20 MHz pumping frequency, the simulations predict an 
EA-limited total error rate of 1.7 x 10^^ for a sinusoidal 
gate drive, and 1.8x10"^ for a square wave drive. Cancel- 
lation of pumping errors of different signs is not included 
in these figures, and hence the relative deviation of the 
average pumped current from the ideal value / = e/ is 
smaller. It remains an open technological challenge to 
realize a similar or better environment for an SET with 
galvanic connections to biasing leads. 

In summary, we have measured single-electron tunnel- 
ing rates in an NIS single-electron box near the degener- 
acy point in two qualitatively different electromagnetic 
environments. The agreement to a model of first or- 
der tunneling activated by absorption of photons from 
a warm environment is excellent. Accounting for the 
enhanced decoupling in the present experiment due to 
lack of galvanic connections, the parameter values we ob- 
tained for the external environment are comparable with 
those presented in Ref. [sj. The results indicate that en- 
vironmental fluctuations deteriorate the performance of 
single-electronic devices due to enhancement of tunnel- 
ing rates in the energetically unfavourable direction. The 
measurements on the sample with an on-chip capacitive 
shunt indicate that the environment-induced error rate of 
the SINIS turnstile can be suppressed to the metrologi- 
cal level with an appropriate filtering of high-frequency 
noise. We emphasize that the EAT model explains the 
observed tunneling rates without any sub-gap states in 
the superconductor. 

I 



for the S and N electrodes, respectively. The average 
energy deposited per individual tunneling event can be 
evaluated as Qf^^s{E)/T{E). For the single-electron box 
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Appendix A: Detector back-action 

Typical value for the average detector current used in 
the real-time counting experiment was 0.5 nA. Low de- 
tector current has the benefit of reduced heating and 
back-action at the cost of signal-to-noise ratio. To as- 
sess experimentally the effect of detector back-action, we 
measured the tunneling rates at base temperature at dif- 
ferent operation points of the detector. In the parameter 
range where tunneling events could still be reliably de- 
tected, we observed that the tunneling rate at the degen- 
eracy point for the shielded sample G2 was modulated 
by factor three. For the unshielded sample Rl, we found 
no change within experimental error. We were able to re- 
produce these results qualitatively using a model similar 
to Ref. 12, where the effective position of the box gate 
fluctuates due to changes in the charge state of the de- 
tector island. Detector back-action is not included in the 
theoretical curves of Figs. 2 and 3, which may affect the 
values of fitted environment parameters for the shielded 
sample, but does not change the conclusions of the paper 
about the role of the external environment. 



Appendix B: Quasiparticle thermalization in the box 
electrodes 

Here, we address the issue of quasiparticle heating 
and thermalization in the electrodes of the hybrid single- 
electron box. The heat associated with quasiparticle tun- 
neling governed by Eq. ([T]) can be calculated as a function 
of the charging energy difference E as 

(Bl) 
(B2) 

I 

with two active charge states, the average heat load is 



Qs{E) = dE, El J^^ dE2 nsiE,)fs{E,)[l - h{E2)]P{Ei - E2 + E), 

Qn{E) ^-^J^^dE, J^^dE2E2nsiE,)fsiEi)[l - ME2)]PiEi - E2 + E). 
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TABLE I. Calculated heat loads to the box electrodes from 
quasiparticle tunneling in various configurations at 60 mK. 



Sample 


NIN state 




NIS state, to N 


NIS state, to S 


Rl 


4.7 X 10"^° 


W 


1.2 X 10"^^ W 


5.3 X 10-2' W 


Gl 


1.7 X 10"^° 


W 


5.5 X 10-2* W 


2.4 X 10-22 W 



given by 



Q 



N,S 



QN,s(^)/r(^) + Qn. si- E)/T{~E) 
1/TiE) + l/n-E) 



(B3) 



Numerical calculations show that the maximum heating 
occurs at the degeneracy point, i.e., E = 0. Using the 
sample and environment parameters extracted from the 
tunneling rate data, the heat loads in samples Rl and Gl 
calculated from Eqs. ()Bip - (|B3|) are presented in Tbl. H 
In the N electrode, hot quasiparticles can thermalize 
through electron-phonon coupling. The electron-phonon 
heat flow is given hy^ 



-ph 



I]V(Ti 



^N,ph) 



(B4) 



where E is the material-dependent electron phonon cou- 
pling constant, V is the island volume, and TN,ph is the 
temperature of the island phonons, assumed to equal 



the sample stage temperature Tq. We use the litera- 
ture value^^ for copper E = 2 x 10^ Wm-'^K^, and 
V = 0.3 X 1.5 X 0.04 fim^. The island phonons are as- 
sumed to be well thermalized to the sample stage. Us- 
ing the heat load values presented in Tbl. U at the base 
temperature Tq = 60 mK, wc find that the N electrode 
quasiparticle temperature is elevated at most 0.02 mK. 

In the S electrode, the electron-phonon heat flow is ex- 
ponentially suppressed below the critical temperature^^. 
Hence, the dominant heat conduction channel is through 
the 0.3 X 1.3 firn^ overlap with the copper shadow. Us- 
ing the specific resistance value 3 kJl /im^ obtained for 
the SET junctions, we get Rt = 8 kJl for the overlap 
junction. We obtain a lower bound for the heat flow by 
considering only first-order quasiparticle tunneling. In 
order to determine the quasiparticle saturation tempera- 
ture, henceforth denoted by Tg , we then evaluate the heat 
flow though the junction from Eq. (|Bip as a function of 
Tg with Tn = 60 mK. Environmental fluctuations can be 
neglected due to large junction capacitance, i. e., we use 
P{E) — S{E). For the S electrode heat loads presented 
in Tbl. U we find that = 116 mK for sample Rl, and 
100 mK for sample Gl. However, the tunneling rates 
in the parameter range accessed in the experiment have 
a very weak dependence on Tg. The theoretical curves 
presented in the article are calculated with Tg = Tq, and 
there is a negligible difference between them and curves 
calculated with Ts = T°. 
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